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Abstract 

In the framework of potential models for heavy quarkonium the mass 
spectrum for the system (be) is considered. Spin-dependent splittings, tak- 
ing into account a change of a constant for effective Coulomb interaction 
between the quarks, and widths of radiative transitions between the (be) 
levels are calculated. In the framework of QCD sum rules, masses of the 
lightest vector B* and pseudoscalar B c states are estimated, scaling rela- 
tion for leptonic constants of heavy quarkonia is derived, and the leptonic 
constant Jb c is evaluated. 

Introduction 

Recently, theoretical interest has been rising to the study of B c meson, the heavy 
(be) quarkonium with open charm and beauty. This interest is stimulated by 
experimental B c search, being performed at FNAL and LEP detectors. 

From one side, similarly to D and B mesons with the open charm and beauty, 
respectively, B c is long-living particle, decaying due to the weak interaction. From 
another side, B c contains the heavy quarks, only, and, therefore, it can be reliably 
described by the use of the methods, developed for the (cc) charmonium and the 
(bb) bottomonium. 

As for the B c production, exact analytic expressions have been recently derived 
for functions of the heavy quark fragmentation into the heavy quarkonium [p], @, |3| , 
in the scaling limit M 2 / s<l. The functions depend on the ratios of the b- and c- 
quark masses to the B c mass. The normalization of the b — ■> B^*'~c fragmentation 
functions also depends on the leptonic constant value for the B c meson. 

In description of the B c decays, it is important to know its spectroscopic 
characteristics. 
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Some preliminary estimates of the bound state masses of the (be) system have 
been made in papers of refs. [|| [5] , devoted to the description of the charmonium 
and bottomonium properties, and in ref.||. Recently in refs.jTJ and ||, revised 
analysis of the B c spectroscopy has been performed in the framework of the 
potential approach and QCD sum rules. 

In the present paper we consider the (be) spectroscopy with taking into the 
account the change of effective Coulomb interaction constant, defining spin- 
dependent splittings of the quarkonium levels. We calculate the widths of ra- 
diative transitions between the levels and analyse the leptonic constant fs c in 
the framework of the QCD sum rules in the scheme, allowing one to derive scaling 
relation for the leptonic constants of the heavy quarkonia. 

In Section 1 we calculate the mass spectrum of the (be) system with the ac- 
count of the spin- dependent forces. In Section 2 the widths of the radiative tran- 
sitions in the B c meson family are evaluated. In Section 3 the leptonic constant 
of B c is calculated. In Conclusion we discuss the obtained results. 

1 Mass Spectrum of B c Mesons. 

The B c meson is the heavy (6c) quarkonium with the open charm and beauty. 
It occupies intermediate place in the mass spectrum of the heavy quarkonia be- 
tween the (cc) charmonium and the (bb) bottomonium. The approaches, applied 
to study the charmonium and the bottomonium, can be expanded to the descrip- 
tion of the B c meson properties, and experimental observation of B c could test 
these approaches and it could be used for the detailed quantitative study of the 
mechanisms of the heavy quark production, hadronization and decays. 

In the present section we obtain the results on the B c meson spectroscopy. We 
will show that below the threshold for the hadronic decay of the (be) system into 
the BD meson pair, there are 16 narrow bound states, decaying, by cascade way, 
into the lightest pseudoscalar B^(0~) state with the mass m(0 _ ) ~ 6.25 GeV. 

1.1 Potential. 

The mass spectra of the charmonium and the bottomonium are studied ex- 
perimentally in details || and they are very well described in the framework 
of the phenomenological potential models of the nonrelativistic heavy quarks 
[|], |5|, [II], [IT], [L3| . To describe the mass spectrum of the (be) system, one would 
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prefer to use the potentials, whose parameters do not depend on the flavours 
of the heavy quarks, composing the heavy quarkonium, i.e. one would use the 
potentials, which rather accurately describe the mass spectra of (cc) as well as 
(bb), with one and the same set of the potential parameters. The use of such 
potentials allows one to avoid an interpolation of the potential parameters from 
the values, fixed by the experimental data on the (cc) and (bb) systems, to the 
values in the intermediate region of the (be) system. 

As it has been shown in ref. [|14[, with accuracy up to an additive shift, the 



potentials, independent of the heavy quark flavours [§], [5], [TO], |TT], |l3l , coincide 
each to other in the region of the average distances between the heavy quarks in 
the (cc) and (bb) systems, so 

0.1 fm < r < 1 fm , (1) 

although those potentials have different asymptotic behaviour in the regions of 
very low (r — > 0) and very large (r — > oo) distances. 

In Cornel model M in accordance with the asymptotic freedom in QCD, the 
potential does behave Coulomb- likely at low distances, and the term, confining 
the quarks, rises linearly at large distances 



Vc(r) = -^ + - 2 +c , (2) 
or or 



so that 



as = 0.36 



a 



2.34 GeV- 1 , 
m c = 1.84 GeV , 

c = -0.25 GeV . (3) 
The Richardson potential [|10]] and its modifications in refs.[|l3j and also 



respects to the behaviour, expected in the framework of QCD, so 

d 3 q ^ 4 12tt 

(2tt) 3 6 3 27 g 2 ln(l + g 2 /A 2 ) 



so that 

A = 0.398 GeV. (5) 
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In the region of the average distances between the heavy quarks (|1|), the 
QCD-motivated potentials allow the approximations by the power (Martin) or 
logarithmic potential. 



The Martin potential has the form [[Tip 

V M {r) = -c M + d M {K M rY 

so that 





= 1 GeV , 


k 


= 0.1, 


m b 


= 5.174 GeV , 


m c 


= 1.8 GeV , 


cm 


= 8.064 GeV , 


d M 


= 6.869 GeV . 



The logarithmic potential is equal to [12 

V L {r) =c L + d L \n(A L r) 

so that 



Al 

m b 
m c 

cl 
d L 



1 GeV, 
4.906 GeV , 
1.5 GeV, 
-0.6635 GW 
0.733 GeV". 



(6) 



(7) 



(8) 



(9) 



The approximations of the nonrelativistic potential of the heavy quarks in the 
distance region (JI|) in the form of the power ([]) and logarithmic (|S]) laws, allow 
one to study its scaling properties. 

In accordance with the virial theorem, the average kinetic energy of the quarks 
in the bound state is determined by the following expression 



<T >- 



1 rdV 

- < —— > 

2 dr 



(10) 
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Then, the logarithmic potential allows one to conclude, that for the quarkonium 

states one gets 

< T L >= const. , (11) 

independently of the flavours of the heavy quarks, composing the heavy quarko- 
nium, 

const. = d L /2^ 0.367 GeV . 

In the Martin potential, the virial theorem ( |I0|) allows one to obtain the 
expression 

< Tm >= Y^k (cm + E) ' (12) 

where E is the binding energy of the quarks in the heavy quarkonium. Phe- 
nomenologically, one has \E\ <C Cm (for example, E(lS,cc) ~ —0.5 GeV), so 
that, neglecting the binding energy of the heavy quarks inside the heavy quarko- 
nium, one can conclude that the average kinetic energy of the heavy quarks is a 
constant value, independent of the quark flavours and the number of the radial 
or orbital excitation. The accuracy of such approximation for < T > is about 
10%, i.e. \AT/T\ ~ 30 40 MeV. 

From the Feynman-Hellmann theorem for the system with the reduced mass 
ft, one has 

£ = -^, (13, 

and, in accordance with condition (|TTD, it follows that the difference of the energies 
for the radial excitations of the heavy quarkonium levels does not depend on the 
reduced mass of the QQ' system 

E(n, n) - E{n, /i) = E{n, //) - E(n, fx') . (14) 

Thus, in the approximation of both the low value for the binding energy of the 
quarks and zero value for the spin-dependent splittings of the levels, the heavy 

Table 1: The Mass Difference for the Two Lightest Vector States of the Different 
Heavy Systems, AM = M (25) - M(1S) in MeV. 



system 


T 




B c 




AM 


563 


588 


585 


660 
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quarkonium state density does not depend on the heavy quark flavours 



dn 

——— = const. (15) 
dM„ 



1 n 



The given statement has been also derived in ref.||16| by the use of the Bohr- 
Sommerfeld quantization of the S-wave states for the heavy quarkonium system 



with the Martin potential [11 



Relations (P^)-([T3|) are phenomenologically confirmed for the vector S-levels 
of the bb, cc, ss systems (see table p]). 

Thus, the structure of the nonsplitted S-levels of the (be) system must not 
only qualitatively, but quantitatively repeat the structure of the S-levels for the 
bb and cc systems, with the accuracy by overall additive shift of the masses. 

Moreover, in the framework of the QCD sum rules, the universality of the 
heavy quark nonrelativistic potential (the independence on the flavours and the 
scaling properties ([TT]) , (p!4|), ([15])) allows one to obtain the scaling relation for 
the leptonic constants of the S-wave quarkonia [|l(| 



f 2 

— = const. , (16) 
independently of the heavy quark flavours in the regime, when 
\m Q - m Q i\ is restricted , Aqqd/ m Q,Q' < 1 , 

i.e., when one can neglect the heavy quark mass difference, and, in the regime, 
when the mass difference is not low, one has 



p /My 

M U/J 



const. , (17) 



where 

_ m Q m QI 
^ m Q + m Q > 



Consider the mass spectrum of the (be) system with the Martin potential |TT 
Solving the Schrodinger equation with the potential (j^) and the parameters 
(0), one finds the B c mass spectrum and the characteristics of the radial wave 
functions -R(O) and R'(0), shown in the tables and [3|, respectively. 
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Table 2: The Energy Levels of the be System, Calculated with no Taking into the 
Account Relativistic Corrections, in GeV. 



n 


1 


0] 




PI 


n 


11 







|15| 


n 


1 


|17| 


B 


IS 


6.301 


6.315 


6.344 


2P 


6.728 


6.735 


6.763 


3D 


7.008 


7.145 


7.030 


2S 


6.893 


7.009 


6.910 


3P 


7.122 






7.160 


4D 


7.308 




7.365 


3S 


7.237 




7.024 


4P 


7.395 








5D 


7.532 







The average kinetic energy of the levels, lying below the threshold for the (be) 
system decay into the BD pair, is presented in table |], from which one can see 
that the term, added to the radial potential due to the orbital rotation, 

f 2 

w 1 



2/ir 2 



weakly influences on the value of the average kinetic energy, and the binding 
energy for the levels with L ^ is essentially determined by the orbital rotation 
energy, that does not still strongly depend on the quark flavours (see table |), so 
that the structure of the nonsplitted levels of the (be) system with L / must 
quantitatively repeat the structure of the charmonium and bottomonium levels, 
too. 



1.2 Spin-dependent Splitting of the (be) Quarkonium. 



In accordance with the results of refs.JI 

orbital and spin-spin interactions, causing the splitting of the nL-levels 



1S|| , to take into the account the spin- 

n is 



Table 3: The Characteristics of the Radial Wave Functions R n s(fy (hi GeV 3 / 2 ) 
and R' nP (0) (in GeV 5 / 2 ), Obtained from the Schrodinger Equation. 



n 


Martin 


i 


Ris(o) 


1.31 


1.28 


R2S(0) 


0.97 


0.99 


R' 2P (0) 


0.55 


0.45 


^p(O) 


0.57 


0.51 
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Table 4: The Average Kinetic and Orbital Energies of the Quark Motion in the 
(be) System, in GeV. 



nL 


IS 


2S 


2P 


3P 


3D 


< T > 


0.35 
0.00 


0.38 
0.00 


0.37 
0.22 


0.39 
0.14 


0.39 
0.29 



the principal quantum number, L is the orbital momentum), one introduces the 
additional term to the potential, so it has the form 

(L-S c L-S b \ ( dV(r) 8 1\ 
— * — * 

+ \ a s — ^- ^ + ^a s S c ■ S b 4tt 5(r) (19) 

3 m c mb r 6 3 3m c mb 

+ ^ "5 (3(5 C • n) (Sj, • n) - S c ■ S b ) \ , ft = - , 

where V(r) is the phenomenological potential, confining the quarks, the first 
term takes into the account the relativistic corrections to the potential V(r); the 
second, third and fourth terms are the relativistic corrections, coming from the 
account of the one gluon exchange between the b and c quarks; as is the effective 
constant of the quark-gluon interaction inside the (be) system. 

The value of the as parameter can be determined by the following way. 

The splitting of the S-wave heavy quarkonium (Q1Q2) is determined by the 
expression 

AM(nS) = a s — ^— \R n ,s(0)\ 2 , (20) 
\)m\m2 



Table 5: The Average Energy of the Orbital Motion in the Heavy Quarkonia, in 
the Model with the Martin Potential, in GeV. 



system 


cc 


be 


bb 


m(2P) 


0.23 


0.22 


0.21 
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where R n s(0) is the value of the radial wave function of the quarkonium, at the 
origin. Using the experimental value of the S-state splitting in the cc system || 

AM(1S, cc) = 117 ±2 MeV , (21) 

and the -Ris(O) value, calculated in the potential model for the cc system, one 
gets the mo del- dependent value of the asfy) constant for the effective Coulomb 
interaction of the heavy quarks (in the Martin potential, one has ots{ip) = 0.44). 

In ref.0 the effective constant value, fixed by the described way, has been 
applied to the description of not only the cc system, but also the be and bb 
quarkonia. 

In the present paper we take into the account the variation of the effective 
Coulomb interaction constant versus the reduced mass of the system (//) . 

In the one-loop approximation at the momentum scale p 2 , the "running" cou- 
pling constant in QCD is determined by the expression 

where b = 11 — 2n//3, and n/ = 3, when one takes into the account the contri- 
bution by the virtual light quarks, p 2 < m 2 



c,b- 



In the model with the Martin potential, for the kinetic energy of the quarks 
(cc) inside xp, one has 

< T ls (cc) >~ 0.357 GeV , (23) 
so that, using the expression for the kinetic energy, 



one gets 



<T>= , (24) 



61n(2 <T> n/Afau) 



so that ats(ip) = 0.44 at 

A QCD ^ 164 MeV . (26) 

As it has been noted in the previous section, the value of the kinetic energy 
of the quark motion weakly depends on the heavy quark flavours, and it, practi- 
cally, is constant, and, hence, the change of the effective as coupling is basically 
determined by the variation of the reduced mass of the heavy quarkonium. In 
accordance with eqs. (p5|)-([26|) and table |], for the (be) system one has 
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nL IS 2S 2P 3P 3D 
a s 0.394 0.385 0.387 0.382 0.383 



Note, the Martin potential leads to the Ris(0) values, which, with the accu- 
racy up to 15 -T- 20%, agrees with the experimental values of the leptonic decay 
constants for the heavy cc and bb quarkonia. The leptonic constants are deter- 
mined by the expression 



4vr 9 , /J 



T(QQ -> t+l-) = -e 2 Q a\ m jf- , (27) 



where cq is the heavy quark charge. 
In the nonrelativistic model one has 



\ 



3 



Ris(0) . (28) 



For the effective Coulomb interaction of the heavy quarks in the basic lS-state 
one has 

/A \3/2 

RC™\0) = 2 n a s ) . (29) 
One can see from table ^ that, taking into the account the variation of the effec- 



tive as constant versus the reduced mass of the heavy quarkonium (see eq. (|25|) ), 
the Coulomb wave functions give the values of the leptonic constants for the 
heavy lS-quarkonia, so that in the framework of the accuracy of the potential 
models, those values agree with the experimental values and the values, obtained 
by the solution of the Schrodinger equation with the given potential. 

The taking into the account the variation of the effective Coulomb interaction 
constant becomes especially notable for the T-particles, for which as(T) ~ 0.33 
instead the fixed value as = 0.44. 

Thus, calculating the splitting of the (be) levels, we take into the account the 
as dependence on the reduced mass of the heavy quarkonium. 

As one can see from eq. (|I~9"D , in contrast to the LS-coupling in the (cc) and 
(bb) systems, there is the jj-coupling in the heavy quarkonium, where the heavy 
quarks have different masses (here, LS C is diagonalized at the given J c momentum, 
(J c = L + S c , J = J c + Sb), J is the total spin of the system). We use the following 
spectroscopic notations for the splitted levels of the (be) system, - n 2jc Lj. 
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Table 6: The Leptonic Decay Constants of the Heavy Quarkonia, the Values, 
Measured Experimentally and Obtained in the Model with the Martin Poten- 
tial, in the Model with the Effective Coulomb Interaction and from the Scaling 
Relation (*), in GeV. 



model 


exp.i 


Martin 


Coulomb 


* 


ftp 


410 ± 15 


547 ± 80 


426 ± 60 


410 ±40 


fs c 




510 ±80 


456 ± 70 


460 ± 60 


fr 


715 ± 15 


660 ± 90 


772 ± 120 


715 ± 70 



One can easily show, that independently of the total spin J projection one 

has 

\ 2L+1 L L+1 > = |J = L + 1, S = 1 > , 

\ 2L ~ 1 L L - 1 > = \J = L-1, S=l>, (30) 



\ 2L+1 L L > = J—^—\J = L,S = l>+J^±\J = L,S = 0> 



1+1 I T T C 1 L 

2TTI |J = L ' S = 1> -V2TTi' 



\ 2L ~ l L L > = \l Z ' ~ \J = L, S = 1 > — \ LS A J = L, S = 0> 



where | J, S > are the state vectors with the given values of the total quark spin 
S = S c + Sb, so that the potential terms of the order of l/m c mb, V m 6 lead, 
generally speaking, to the mixing of the levels with the different J c values at the 
given J values. The tensor forces (the last term in eq.flTQD) are equal to zero at 
L = or S = 0. 

One can easily show, that 

3(n p n q - ^5 pq )S p c S q b = ^(n p n q - -5 pq )S p S q , (31) 
since for the quark spin one has 

SSSJ + SJSJ = l Sm ■ (32) 
The averaging over the angle variables can be represented in the form 

< L, m\n p n q \L, m' >= a(L p L q + L q L p ) mm , + b5 pq , (33) 
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where L are the orbital momentum matrices in the respective irreducible repre- 
sentation. 

Let us use the following conditions. 

1) The normalization of the unit vector, 

< n p n q > 5 pq = 1 . (34) 

2) The orthogonality of the radius-vector to the orbital momentum, 

n p L p = , (35) 

3) The commutation relations for the angle momentum, 

[L p - L q ] = ie pql U . (36) 

Then one can easily find, that in eq. (|3~3"D one gets 

1 



4L 2 - 3 ' 
2L 2 - 1 



4L 2 -3 

Thus, (see also ref.fll9|) 



(37) 
(38) 



< 6(nV - y pq )S p S q b >= - 1 3 (6(Lg) 2 + 3{LS) - 2L 2 S 2 ) . (39) 

Using eqs.(^OD) (|39|), for the level shift, calculated in the perturbation theory 
at S — 1, one gets the following formulae 

A£ n i 5o = -a s \R nS (0)\ 2 , (40) 
3m c m{, 

AE n i Sl = a s —^—\R nS (0) | 2 , (41) 

6 1 

AE n sp 2 = a s < > + 

-(— + —) < + - a — > (42) 

4 \m? c m 2 J rdr 3 r 3 

4 1 

AE n i Po = -a s < — > - 

m c mb r d 
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7.0 



6.0 



3S 


AP 


2+ 


AD 

BD threshold 






1'+ " 








1+ 








0+ 


3D 








i- 










2S 






2" 










0" 








2P 


2+ 






1'+ 








1+ 








0+ 




IS 


1- 








0" 







Figure 1: The Mass Spectrum of the (be) System with the Account of the Split- 
tings. 

1/1 1 \ dV(r) 8 1 

2 m\) rdr 3 S r 3 ' 

52 1 

A£„5 D3 = a s — < -3 > + 

2lm c rrib r 6 

1/1 1 \ dV ( r ) 8 1 

2U + W) < ^ + 3 a V > ' (44) 

92 1 

AE n s Dl = -a s — < -r > - 

21m c rrib r 6 
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3/1 M <W(r) 81 



C 



where _R n s(0) are the radial wave functions at L = 0, < ... > denote the average 
values, calculated under the wave functions i? n ^(r). The mixing matrix elements 
have the forms 

< 3 P 1 \AE\ 3 P 1 > = -a s —^- + 

\)m c rrib r 6 

_1 5_\ < dV(r) + 8 ^ 1_ > ^ 4g ^ 

Ami 12ml J r ^ r 3 r 3 

9m c mi) r 6 



1 + J_\ < dVjr) + 8 ^ 1_ > 
2m^ 6ml J rdr 3 r 3 



, , , , 2^ 1 

< 3 P X A£ ^ > = -e* s — 1— < - > 

\)m c mb r 6 



V2 ^(r) 8 1 

< it 2 + 77 as-, > , (4* 



6m^ rdr 3 r 3 

4 1 

< 5 D 2 \AE\ 5 D 2 > = -a s — < -r > + 



15m c rrib r 3 

_1 1_\ < dVjr) + 8 ^ 1_ > ^ 4g ^ 

2m\ 5ml J rdr 3 r 3 



< 3 D 2 |AP| 3 P 2 > = -a s — — < — > + 

lom c m& r 6 



3 9 \ dV(r) 8 1 

S <-^ + o«5-r>, (50) 



Am 2 c 20ml ' r ^ r 

< 5 D 2 |AP| 3 P 2 > = - as < 1 > 

lbm c mb r 6 



V6 dV{r) 8 1 

< r^ + oas^>, ( 51 ) 



10ml rdr 3 r 3 

As one can see from eq.([|l|), the S-level splitting is essentially determined by 
the |i?ns(0)| value, which can be related with the leptonic decay constants of the 
S-states (CT, 1~). Section 3 is devoted to the calculation of these constants by 
different ways. We note here only, that with enough accuracy, the predictions of 
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different potential models on the |-Ris(0)| value are in an agreement with each 
other as well as with predictions in other approaches. 

For the 2P, 3P and 3D levels, the mixing matrices of the states with the total 
quark spin S = 1 and S = have the forms 

\2P, l'+ > = 0.294|S = 1 > +0.95615" = > , 
\2P, l + > = 0.95615 = 1 > -0.29415" = > , 

so that in the 1 + state the probability of the total quark spin value S 
to w{2P) = 0.913, 

1 > +0.929|5 = > , (54) 
1 > — 0.371|5 = > , (55) 

= 1 > +0.825|S = > , (56) 
1 > +0.566|5 = > , (57) 

so that w(3D) = 0.680. 

With the account of the calculated splittings, the B c mass spectrum is shown 
on fig.|I| and in table [7| 

The masses of the B c mesons have been also calculated in papers of ref. |28| . 

As one can see from tables [2] and [5], the place of the lS-level in the (bc) 
system (m(lS) ~ 6.3 GeV) is predicted by the potential models with the rather 
high accuracy Am(lS) ~ 30 MeV, and the lS-level splitting into the vector and 
pseudoscalar states is about m(l~) — m(0~) ~ 70 MeV. 

1.3 B c Meson Masses from QCD Sum Rules. 

Potential model estimates for the masses of the lightest (be) states are in the 
agreement with the results of the calculations for the vector and pseudoscalar 
(be) states in the framework of the QCD sum rules [§, |29], |3T| , where the calcu- 
lation accuracy is lower, than the accuracy of the potential models, because the 
results essentially depend on both the modelling the nonresonant hadronic part 
of the current correlator (the continuum threshold) and the parameter of the sum 
rule scheme (the number of the moment for the spectral density of the current 
correlator or the Borel transformation parameter), 

m SR (0~) « m SR (l~) ~ 6.3 6.5 GeV . (58) 



(52) 
(53) 

1 is equal 



\3P, 1'+ > = 0.37115 = 

\3P, 1+ > = 0.929|5 = 

so that w(3P) = 0.863, 

\3D, 2'- > = -0.566|5 = 

|3L>, 2" > = 0.825|5 = 
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Table 7: The Masses (in GeV) of the Bound (be) States below the Threshold of 
the Decay into the BD Meson Pair, * is the present paper. 



state 


* 


E 


115 




1% 


6.253 


6.264 


6.314 


1% 


6.317 


6.337 


6.355 


2% 


6.867 


6.856 


6.88 


9 


2 1 Si 


6.902 


6.899 


6.917 


2 l P 


6.683 


6.700 


6.72 


8 


2P 1+ 


6.717 


6.730 


6.760 


2P 1 /+ 


6.729 


6.736 






2 3 P 2 


6.743 


6.747 


6.773 


&P 


7.088 


7.108 


7.134 


3P 1+ 


7.113 


7.135 


7.159 


3P 1'+ 


7.124 


7.142 






3 3 P 2 


7.134 


7.153 


7.166 


3P 2~ 


7.001 


7.009 






3 5 P 3 


7.007 


7.005 






3 3 Pi 


7.008 


7.012 






3P 2 /_ 


7.016 


7.012 







As it has been shown in papers of ref.||31||, for the lightest vector quarkonium, 
the following QCD sum rules take place 

_ftM$_ = l p ^mIif D ^\s) ^ + u Q CDnonpert{q2) ^ (5g) 

TTly — Q VT J Si S — Q 



where fy is the leptonic constant of the vector (be) state with the mass M< 



v 



zf v M v e x / px = <0\J,(x)\V(p, A)>, (60) 
J^(x) = c(a:)7 M 6(x) , (61) 

where A, p are the B* polarization and momentum, respectively, and 

Jd 4 xe^ <0\TJ,(x)J u (0)\0> = (-g^ + ^j n$ CD + qil q u Tl Q s CD (Q2) 
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Table 8: The Masses (in GeV) of the Lightest Pseudoscalar B c and Vector B* 
States in Different Models, * is the present paper. 



state 


* 


|17| 


[0 





11 


|21| 


H E 


o- 

1" 


6.253 
6.317 


6.249 
6.339 


6.314 
6.354 


6.293 
6.346 


6.270 
6.340 


6.243 
6.320 


6.246 
6.319 


state 


11 


11 


m 


HI 


H 


|2§ 


§ 


o- 

1- 


6.264 
6.337 


6.320 
6.370 


6.256 
6.329 


6.276 
6.365 


6.286 
6.328 


6.320 


6.255 
6.330 



TlQ CD (q 2 ) = Ttf CDpert + nQCDnonpert^ ? (63) 
u QCDnon P ert^ = ^C^O', (64) 

where O l are the vacuum expectation values of the composite operators such as 
< rmpip >, < as G 2 > and so on. The Wilson coefficients are calculable in 
the perturbation theory of QCD. = (m c + nib) 2 is the kinematical threshold 
of the perturbative contribution, M v > Sj, s t h is the threshold of the nonreso- 
nant hadronic contribution, which is considered to be equal to the perturbative 
contribution at s > s t h- 

Considering the respective correlators, one can write down the sum rules, 
analogous to eq. ([59"D , for the scalar and pseudoscalar states. 

One believes that the sum rule ( |59D must rather accurately be valid at q 2 < 0. 

For the n-th derivative of eq.(|59"D at q 2 = one gets 

#<«?r = I jT amn !r"' (s) * + ^ ^n?~'(^) , (65) 

so, considering the ratio of the n-th derivative to the n+ 1-th one, one can obtain 
the value of the vector B* meson mass. The result of the calculation depends on 
the n number in the sum rules (|65|) , because of taking into the account both still 
finite number of the terms in the perturbation theory expansion and restricted 
set of the composite operators. 

The analogous procedure can be performed in the sum rule scheme with the 
Borel transform, leading to the dependence of the results on the transformation 
parameter. 
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As one can see from eq.(|B5D, the result, obtained in the framework of the 
QCD sum rules, depends on the choices of the values for the hadronic continuum 
threshold energy and the current masses of the quarks. Then, this dependence 
causes large errors in the estimates of the masses for the lightest pseudoscalar, 
vector and scalar (be) states. 

Thus, the QCD sum rules give the estimates of the quark binding energy in 
the quarkonium, and the estimates are in the agreement with the results of the 
potential models, but they contain the large parametric uncertainty. 



2 Radiative Transitions in the B c Family. 

The B c mesons have no annihilation channels for the decays due to QCD and elec- 
tromagnetic interactions. Therefore, the mesons, lying below the threshold for 
the production of B and D mesons, will, by cascade way, decay into the 0~(1S') 
state by emission of 7 quanta and 7r mesons. Theoretical estimates of the tran- 
sitions between the levels with the emission of the ir mesons have uncertainties, 
and the electromagnetic transitions are rather accurately calculable. 

2.1 Electromagnetic Transitions. 



The formulae for the radiative El-transitions have the form |32|, |33 

4 



T(nPj - 


,n 1 S 1 + 1 ) = 


T(nPj - 




T{n 1 S 1 - 





dem Qlf f oo 3 I 2 (nP; nS) Wj(nP) 

4 



a em Q 2 eff oo 3 I 2 (nP; nS) (1 - Wj(nP)) 

^a em Q 2 eff u 3 I 2 (nS;nP) (2 J + 1) 

wj{nP) , (66) 
4 
9 

1 - Wj(nP)) 
4 

2"7 

wj(nP))wj>(nD)Sjj' 
4 



T(n 1 S ^nPj + 1 ) = -a em Q 2 eff u; 3 I 2 (nS;nP) (2J+1) 



V{nPj -> nDj, + 7) = —a em Q 2 eff uj 3 I 2 (nD-nP)(2J' + l 



T(nDj -> nPj, + 7) = ^ a em Q 2 eff u 3 I 2 {nD; nP) (2J' + 1) 
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w J/ (nP))wj(nD)Sj>j , 



where u is the photon energy, a em is the electromagnetic fine structure con- 
stant, wj(nL) is the probability that the spin S = 1 in the nL state, so that 
w (nP) = w 2 (nP) = 1, w\(nD) = w 3 (nD) = 1, and the wi(nP), w 2 {nD) values 
are presented in the previous section. 



The statistical factor Sjj> takes the values [ 33 1 



J J' S 




1 
1 

2 
2 
2 



1 
1 

2 
1 

2 
3 



jj> 
2 

1/2 

9/10 

1/50 

9/50 

18/25 



The I(nL\ nL') value is expressed through the radial wave functions, 



I(nL;nL') = \ J R n L{r)RnL> 



r)r 3 dr\ 



(67) 



For the set of the transitions one obtains 



/(IS, 2P) 
/(IS, 3P) 
/(2S, 2P) 
/(2S, 3P) 
/(3P,2P) 
/(3P,3P) 



1.568 GeV' 1 
0.255 GeV' 1 
2.019 GeV" 1 
2.704 GeV' 1 
2.536 GeV" 1 
2.416 GeV- 1 



(6£ 



In eq. 



one uses 



Q eff = (m c Qi - m b Q c )/(m c + m b ) 



(69) 



where Q c ^ are the electric charges of the quarks. For the B c meson with the 
parameters from the Martin potential, one gets Q e ff — 0.41. 
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For the dipole magnetic transitions one has ||, [32|, |33| 

1 fi 

TtfSi - n'Sf + 7) = - £ ft u; 3 (2/ + 1) A% , (70) 



where 
and 



A if = J R nS (r)Rns(r)jo(ujr/2)r 2 dr 



Veff = O ^ cl7lb ~ Qb^c) ■ (71) 

Note, in contrast to the and T-particles, the total width of the B* meson is 
equal to the width of its radiative decay into the £> c (0~) state. 

The electromagnetic widths, calculated under eqs.(|66D,([70|), and the frequen- 
cies of the emitted photons are presented in tables [LOj, |ll|. 

Thus, the registration of the cascade electromagnetic transitions in the (be) 
family can be used for the observation of the higher (fec)excitations, having no 
annihilation channels of the decays. 



2.2 Hadronic Transitions. 

In the framework of QCD the consideration of the hadronic transitions between 
the states of the heavy quarkonium family is built on the basis of the multipole 
expansion for the gluon emission by the heavy nonrelativistic quarks |[34|| , with 
the postcoming hadronization of the gluons, independently of the heavy quark 
motion. 

In the leading approximation over the velocity of the heavy quark motion, the 
action, corresponding to the heavy quark coupling to the external gluon field, 

Sint = -9 J d*x Al(x) ■ f a {x) , (72) 
can be expressed in the form 

S m t = 9 Jdt r k E a k (t, x) ^m n (r)^ f \r) K(s n , f) d 3 f , (73) 

where ^(^O is the wave function of the quarkonium, emitting the gluon, ^fj(f) is 
the wave function of the colour-octet state of the quarkonium, K(s n , f) respects to 
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the spin factor (in the leading approximation, the heavy quark spin is decoupled 
from the interaction with the gluons). 

Then the matrix element for the El-El transition of the quarkonium nLj — > 
n'L'j, + gg can be written in the form 

M(nLj^ n'L'j, +gg) = 4vra 5 E a k E b m ■ 

J d 3 rd 3 r' r k r' m Gf^r') * nLj {r)^ n , L>ji (r')(74) 

where G^ b _ Sn (r, r') respects to the propagator of the colour-octet state of the heavy 
quarkonium 

G = 73^. < 75 > 

where Hqq is the hamiltonian of the coloured state. 



One can see from eq. (|74|) , that the determination of the transition matrix 
element depends on both the wave function of the quarkonium and the hamilto- 
nian Hqq. Thus, the theoretical consideration of the hadronic transitions in the 
quarkonium family is model dependent. 

In the set of papers of ref. for the calculation of the values such as ([74]) , 
the potential approach has been developed. 



In papers of ref.|SB|] it has been shown that nonperturbative conversion of 
the gluons into the n meson pair allows one to make the consideration in the 
framework of the low-energy theorems in QCD, so that this consideration agrees 
with the papers, performed in the framework of PCAC and soft pion technique 

However, as it follows from eq.(74]) and the Wigner-Eckart theorem, the dif- 



ferential width for the El-El transition allows the representation in the form 
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dm 2 




<11 nLj^n'L' J , + h) = (2J l + l)Y / { ™ " % \ A k (L,L') , (76) 



where m 2 is the invariant mass of the light hadron system h, { } are 6j-symbols, 
A k (L,L') is the contribution by the irreducible tensor of the rang, equal to k = 
0, 1, 2, s is the total quark spin inside the quarkonium. 
In the limit of soft pions, one has A±(L, V) = 0. 
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From eqs.(T^), ( |7"6|) it follows, that, with the accuracy by the difference in 
the phase spaces, the widths of the hadronic transitions in the (QQ) and (QQ') 
quarkonia are related by the following expression [R3, BH| 



r(QQ') < r 2 (QQ') > 5 



(77) 



T(QQ) < r 2 (QQ) > 2 ' 
Then the experimental data on the transitions of ip' — > J /ip + tttt, T' — > T + 7r7r 



■^(3770) — > J/^> + 7T7T |38j allow one to extract the values of Ak(L,L') for the 
transitions 2S* — > IS + 7T7t and 3D — > IS* + 7T7t [[?]]. 

The invariant mass spectrum of the 7r meson pair has the universal form 



36, 37 



where a; = m/2m 7r , | &„.„. | is the tttt pair momentum. 

The estimates for the widths of the hadronic transitions in the (be) family 
have been made in ref.0. The hadronic transition widths, having the values 
comparable with the electromagnetic transition width values, are presented in 
table [l^. The transitions in the (be) family with the emission of r] mesons are 
suppressed by the low value of the phase space. 

Thus, the registration of the hadronic transitions in the (be) family with the 
emission of the it meson pairs can be used to observe the higher 2S- and 3D- 
excitations of the basic state. 



3 Leptonic Constant of B c Meson. 

As we have seen in Section 1, the value of the leptonic constant of the B c meson 
determines the splitting of the basic lS-state of the (be) system. Moreover, the 
higher excitations in the (be) system transform, by cascade way, into the lightest 
0~ state of B c , whose widths of the decays are essentially determined by the 



value of fs c , too. In the quark models [[39], £40], [|T|, used to calculate the weak 
decay widths of mesons, the leptonic constant, as the parameter, determines the 
broadness of the quark wave packet inside the meson (generally, the wave function 
is chosen in the oscillator form), therefore, the practical problem for the extraction 
of the value for the weak charged current mixing matrix element \Vb c \ from the 
data on the weak B c decays can be only solved at the known value of fs c - 
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Thus, the leptonic constant fs c is the most important quantity, characterizing 
the bound state of the (be) system. 

In the present Section we calculate the value of fg c in different ways. 

To describe the bound states of the quarks requires the use of the nonpertur- 
bative approaches. The bound states of the heavy quarks allow one to consider 
simplifications, connected to both large values of the quark masses ^qcd/i^q *C 1 
and the nonrelativistic quark motion v — *■ 0. Therefore the value of fs c can be 
rather reliably determined in the framework of the potential models and the QCD 
sum rules |3l| . 



3.1 Jb c from Potential Models. 

In the framework of the nonrelativistic potential models, the leptonic constants 
of the pseudoscalar and vector mesons (see eqs.(^), (j6T|) ) 

<0\c(x)^b(x)\B* c (p,e) > = if v M v e ipx , (79) 
<0\c(x) l5lfl b(x)\B c (p) > = z/pjve^, (80) 



are determined by expression (|28] 

fv = fp 



\ 7iM Bc{ls) 



Ris(o) 



where Ris{0) is the radial wave function of the lS-state of the (be) system, at the 
origin. The wave function is calculated by solving the Schrodinger equation with 
the different potentials [f|, |5], [03' O' HH H3' ^ n ^ ne quasipotential approach |42] or 



by solving the Bethe-Salpeter equation with instant potential and in expansion 
up to second order over the quark motion velocity v/c ||3|, pHfl . 

The values of the leptonic B c meson constant, calculated in the different po- 
tential models and effective Coulomb potential with the "running" as constant, 
determined in Section 1, are presented in table [L3[ 

Thus, in the approach accuracy, the potential quark models give the fs c 
values, being in the good agreement with each other, so that 

fP°* = 500 ± 80 MeV . (82) 
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3.2 f Bc from QCD Sum Rules. 



In the framework of the QCD sum rules |51], expressions d59|)-(|65|) have been 
derived for the vector states. The expressions are considered at q 2 < in the 
schemes of the spectral density moments QS5| ) or with the application of the Borel 
transform f3~If . As one can see from eqs.(|59D - (|65|), the result of the QCD sum rule 
calculations is determined by not only the physical parameters such as the quark 
and meson masses, but also by the unphysical parameters of the sum rule scheme 
such as the number of the spectral density moment or the Borel transformation 
parameter. In the QCD sum rules, this unphysical dependence of the fs c value is 
caused by that the consideration is yet performed with the finite number of terms 
in the expansion of the QCD perturbation theory for the Wilson coefficients of 
the unit and composite operators. In the calculations, the set of the composite 
operators is also restricted. 

Thus, the ambiguity in the choice of the hadronic continuum threshold and the 
parameter of the sum rule scheme essentially reduces the reliability of the QCD 
sum rule predictions for the leptonic constants of the vector and pseudoscalar B c 
states. 

Moreover, the nonrelativistic quark motion inside the heavy quarkonium 
v — > leads to that the a^/f-corrections to the perturbative part of the 
quark current correlators become the most important, where as is the effective 
Coulomb coupling constant in the heavy quarkonium. As it has been noted in 
refs.[[U], [H], £|6|], the Coulomb as/v -corrections can be summed and represented 
in the form of the factor, corresponding to the Coulomb wave function of the 
heavy quarks, so that 

F{v) = ^ ^ , (83) 

v ; 3v 1 -exp(-47ra s /3v) v ; 

where 2v is the relative velocity of the heavy quarks inside the quarkonium. The 
expansion of the factor fl83f) in the first order over as/v 



F(v) ~ 1 - ^p- , (84) 
Sv 

gives the expression, obtained in the first order of the QCD perturbation theory 



Note, the as parameter in eq.(p3|) has to be at the scale of the characteristic 
quark virtualities in the quarkonium (see Section 1), but at the scale of the quark 
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or quarkonium masses, as sometime one does it and decreases the value of the 
factor (|83|) . 

The choice of the as parameter essentially determines the spread of the sum 
rule predictions for the /b c value (see table [14]) 

ff* = 160 -T- 570 MeV . (85) 



As one can see from eq.(f^), the ambiguity in the choices of the QCD sum rule 
parameters leads to the essential deviations of the results from the /g c estimates 
( |8~2"D in the potential models. 

However, as it has been noted in Section 1, 

i ) the large value of the heavy quark masses Aqcd/^iq "C 1, 

ii ) the nonrelativistic heavy quark motion inside the heavy quarkonium v — *■ 0, 

and 

iii) the universal scaling properties of the potential in the heavy quarkonium, 
when the kinetic energy of the quarks and the quarkonium state density do 
not depend on the heavy quark flavours (see eqs.( |IO|) - fllB])), 

allow one to state the scaling relation ( |17|) for the leptonic constants of the S-wave 
quarkonia 

p /My 



M V4/i/ 

Indeed, 



const. 



i ) at KqcD/ m Q <^ 1 one can neglect the quark-gluon condensate contribution, 

having the order of magnitude 0(l/m&m c ) (the contribution into the if) and 
T leptonic constants is less than 15%), 

ii ) at v — > one has to take into the account the Coulomb-like a^/f-corrections 

in the form of the factor flB3f) , so that the imaginary part of the correlators 
for the vector and axial quark currents has the form 

QmU v {q 2 ) ~ QmU P {q 2 ) = ^ q 2 ( ^-Y , (86) 



4/i/ 



where 



2 t Ambm c 



q 2 — (nib — trie) 2 
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Moreover, condition ( |T5| ) can be used in the specific QCD sum rule scheme, so 
that this scheme excludes the dependence of the results on the parameters such 
as the number of the spectral density moment or the Borel parameter. 

Indeed, for example, the resonance contribution into the hadronic part of the 
vector current correlator, having the form 



4™ V) = / 7^3 E fkMvAs - My n ) , (87) 

J S Q , r , 



can be rewritten as 



^, (f)= j^ sfv< ^±.^ein-k). (88) 

where n(s) is the number of the vector S-state versus the mass, so that 

n(m 2 k ) = k . (89) 
Taking the average value for the derivative of the step-like function, one gets 

nt-'tf)=<|i:»(»-*)>/ J ^^,^, (") 

and, supposing 

< j-E^( n - k ) >- 1 > ( 91 ) 

k 

one can, in average, write down 

3m < U {hadr \q 2 ) >= 3mn (QCT) (? 2 ) , (92) 

so, taking into the account the Coulomb factor and neglecting power corrections 
over 1/rriQ, at the physical points s n = M% one obtains 



1L f—Y - dM - 



(93) 



M n \4jiJ n dn 
where one has supposed that 

m b + m c ^ M Bc , (94) 
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and 

fvn ^ fpn = fn ■ (95) 

Further, as it has been shown in Section 1, in the heavy quarkonium the value 
of dn/dM n does not depend on the quark masses (see eq. ([l5|) ), and, with the 
accuracy up to the logarithmic corrections, a$ is the constant value (the last 
fact is also manifested in the flavour independence of the Coulomb part of the 
potential in Cornel model). Therefore, one can draw the conclusion that, in the 
leading approximation, the right hand side of eq.(p3|) is the constant value, and 
there is the scaling relation ( |T7D [[LlJ . This relation is valid in the resonant region, 
where one can neglect the contribution by the hadronic continuum. 

Note, scaling relation ( |T7j ) is in the good agreement with the experimental 
data on the leptonic decay constants of the ip- and T-particles (see table ^), for 
which one has 4/z/M = 1 [|TJ| . 

The value of the constant in the right hand side of eq. fllTD is in the agreement 
with the estimate, when we suppose 

dMr 1 

< >~ -((Mr - M r ) + (M T » - Mr')) , (96) 

dn 2 

and as = 0.36, as it is in Cornel model. 

Further, in the limit case of B- and .D-mesons, when the heavy quark mass is 
much greater than the light quark mass niQ ^> m q , one has 

u ~ m q 

and 

7r dn 

Then it is evident that at one and the same u one gets 

f M = const. (98) 



Scaling law ([5B|) is very well known in EHQT |St| for mesons with a single heavy 
quark (Qq), and it follows, for example, from the identity of the B- and D-meson 
wave functions in the limit, when infinitely heavy quark can be considered as a 
static source of gluon field. 

In our derivation of eqs. (|97|) and ( |9~8|) we have neglected power corrections over 
the inverse heavy quark mass. Moreover, we have used the presentation about 
the light constituent quark with the mass, equal to 

m q ~ 330 MeV , (99) 
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so that this quark has to be considered as nonrelativistic one v — > 0, and the 
following conditions take place 

mg + m,~ M^ q) , m q < m Q , (100) 

and 

fv^fp = f- (101) 
In agreement with eqs. (|97D and (0), one finds the estimates^ 

f BM = 120 ± 20 MeV , (102) 
/dm = 220 ± 30 MeK , (103) 

that is in an agreement with the estimates in the other schemes of the QCD sum 



rules [31, 51 



Thus, in the limits of 4/x/M = 1 and fi/M -C 1, scaling relation (|17]) is 
consistent. 

The fB c estimate from eq. ([l7D contains the uncertainty, connected to the 
choice of the ratio for the b- and c-quark masses, so that (see table |b|) 

f Bc = 460 ± 60 MeV . (104) 

In ref. |^6[ the sum rule scheme with the double Borel transform has been used. 
So, it allows one to study effects, related to the power corrections from the gluon 
condensate, corrections due to nonzero quark velocity and nonzero binding energy 
of the quarks in the quarkonium. 

Indeed, for the set of narrow pseudoscalar states, one has the sum rules 

I (mt+m M % g -^) = y^f amiws) + ca(i) < > • < io5 > 

where 

GKq J 192m b m c q 2 \ 2v 5 1-v v 4 /' V ' 



1 In ref. [^6| the dependence of the S-wave state density dn/dM n on the reduced mass of the 
system with the Martin potential has been found by the Bohr-Sommerfeld quantization, so that 
at the step from (bb) to (bq), the density changes less than about 15%. 



28 



and 

f = q 2 - (m b - m c y , v 2 = 1 _W^c_ (1Q7) 

Acting by the Borel operator L T (—q 2 ) on eq. ( |105|) , one gets 
r4f2 

kJ PI 

^ (m& + m c y 
where 



y m kIPk e -ulr = L\ ds ^ m U P (s) e~ ST + C' g (t) < —G 2 > , (108) 
r-| ra h + m c r vr J tt 



L T (x) = lim — - (-— , n/x = r, (109) 



n,a;-»oo n \ \ dx 
C' g {t) = L T (-q 2 ) C G (q 2 ) . (110) 

For the exponential set in the left hand side of eq.( |108|) , one uses the Euler- 
MacLaurin formula 



E M k JPk „-MIt f°° J-A/r ,,4,2 -Mir , \ ~> » /f 4 f 2 „-M%Tj_ 
— e k = dM k ——M k f Pk e * + }^ M kJpk e k +•••■ 

fe=l ( m & + m cr -W dMfc ^ 

(111) 

Making the second Borel transform L M |(r) on eq. ( |108| ) with the account of 
eq. (|lll|) , one finds the expression for the leptonic constants of the pseudoscalar 
(be) states, so that 

fh = 2( "y' iJ w + c 'i^ < v G ' 2 > } • (112) 

where we have used the following property of the Borel operator 

L T {x) x n e~ bx -> ^(r - b) . (113) 

Explicit form for the spectral density and Wilson coefficients can be found in 
ref.ll ■ 

Expression ( |112j ) is in the agreement with the above performed derivation of 
scaling relation fll7|) . 

The numerical effect from the mentioned corrections considers to be not large 
(the power corrections are of the order of 10%), and the uncertainty, connected 
to the choice of the quark masses, dominates in the error of the fs c value deter- 
mination (see eq. (|104j) ). 

Thus, we have shown that, in the framework of the QCD sum rules, the most 
reliable estimate of the fs c value ( |104| ) is coming from the use of the scaling 
relation (|17|) for the leptonic decay constants of the quarkonia, and this relation 
very well agrees with the results of the potential models. 
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Conclusion. 



In the present paper we have considered the spectroscopic characteristics of the 
bound states in the (be) system. 

We have shown that below the threshold of the (be) system decay into the BD 
meson pair, there are 16 narrow states of the B c meson family, whose masses can 
be reliably calculated in the framework of the nonrelativistic potential models of 
the heavy quarkonia. The flavour independence of the QCD-motivated poten- 
tials in the region of average distances between the quarks in the (bb), (cc) and 
(bc) systems and their scaling properties allow one to find the regularity of the 
spectra for the levels, nonsplitted by the spin-dependent forces: in the leading 
approximation the state density of the system does not depend on the heavy 
quark flavours, i.e. the distances between the nL-levels of the heavy quarkonium 
do not depend on the heavy quark flavours. 

We have described the spin-dependent splittings of the (be) system levels, so, 
the splittings, appearing in the second order over the inverse heavy quark masses, 
Vsd — 0(l/m b m c )i with the account of the variation of the effective Coulomb 
coupling constant of the quarks (the interaction is due to relativistic corrections, 
coming from the one gluon exchange). 

The approaches, developed to describe emission by the heavy quarks, have 
been applied to the description of the radiative transitions in the (be) family, 
whose states have no electromagnetic or gluonic channels of the annihilation. 
The last fact means that, due to the cascade processes with the emission of 
photons and pion pairs, the higher excitations take the transitions into the lightest 
pseudoscalar B c meson, decaying by the weak way. Therefore, the excited states 
of the (be) system have the widths, essentially less (by two orders of magnitude) 
than the widths in the charmonium and bottomonium systems. 

As for the value of the leptonic decay constant /b c , it can be the most reliably 
estimated from the scaling relation for the leptonic constants of the heavy quarko- 
nia, due to the relation, obtained in the framework of the QCD sum rules in the 
specific scheme. In the other schemes of the QCD sum rules, it is necessary to 
do an interpolation of the scheme parameters (the hadronic continuum threshold 
and the number of the spectral density moment or the Borel parameter) into the 
region of the (be) system, so this procedure leads to the essential uncertainties. 
The fs c estimate from the scaling relation agrees with the results of the potential 
models, whose accuracy for the leptonic constants is notably lower. The value of 
fs c essentially determines the decay widths and the production cross sections of 
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the B c mesons. 

The B c decays have been studied in refs. 0, |29], [52], where it has been shown, 
that the B c lifetime is approximately equal to 

t(B c ) ~ 0.5^0.7 ps , (114) 

and the characteristic decay mode, having the preferable signature for the exper- 
imental search, has the branching fraction, equal to 

Br(B+ -> ipX) ~ 17% . (115) 

The decay of — > ipir + , having the branching fraction 

Br{B+ -»• ~ 0.2% , (116) 

is chosen by CDF Collaboration for the B c search, when about 20 events are 
expected |53fl . 

The B c production at the FNAL and LEP colliders has been studied in refs. 
0, 0, H |53], |54j], where it has been shown, that the level of the B c yield has the 
order of 

a(B c X) , 

Thus, in the present paper we have described the spectroscopic characteristics 
of the B c mesons, whose search at FNAL and LEP is carrying out, and, as 
expected, it will be successfully realized in the nearest future. 
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Table 9: The Energies (in MeV) and Widths (in keV) of the Electromagnetic 
El-Transitions in the (be) Family. 



transition 
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101 -|- y 
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1 31 1 
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' 100 1 7 


464 


11.6 


0.0 
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1 Q 9 
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3P ^ 
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Table 10: The Energies (in MeV) and Widths (in keV) of the Electromagnetic 
El-Transitions in the (be) Family. 



transition 


CO 


1 


1 [ill 


3P 2 — ► 


3iJi + 7 


IzO 


n 1 
U.l 


U.z 


O/2 — > 


0'— 1 ~, 
01J z +7 


1 1 fi 
1 lo 


u.o 




3P 2 — > 


oiJ z +7 


loo 


1.0 


O.Z 


Q P, > 
O.T2 ' 


Oi/3 -f- 7 


1Z ( 


inn 
iu.y 


1 7 s 
1 1 .0 




3Dx + 7 


80 


3 2 


6 9 


3P 1'+ 


-> 3£>! + 7 


116 


0.3 


0.4 


3P 1+ 


-> 3Di + 7 


105 


1.6 


0.3 


3P 1'+ 


-> 3P> 2'- + 7 


108 


3.5 




3P 1+ 


-> 3P> 2- + 7 


112 


3.9 


9.8 


3P 1'+ 


-> 3P> 2- + 7 


123 


2.5 


11.5 


3P 1+ 


-> 3D 2'- + 7 


97 


1.2 




3D 3 ^ 


2P 2 + 7 


264 


76.9 


98.7 


3P>i -» 


2P + 7 


325 


79.7 


88.6 


3P>i -» 


2P 1'+ + 7 


279 


3.3 


0.0 


3P>i -» 


2P 1+ + 7 


291 


39.2 


49.3 


3P>! -» 


2P 2 + 7 


265 


2.2 


2.7 


3P> 2'- 


^2P 2 + 7 


273 


6.8 




3D 2'- 


-> 2P 2 + 7 


258 


12.2 


24.7 


3D 2'- 


-> 2P 1'+ + 7 


287 


46.0 


92.5 


3D 2'- 


-> 2P 1+ + 7 


301 


25.0 




3D 2- 


-> 2P 1'+ + 7 


272 


18.4 


0.1 


3D 2- 


-> 2P 1+ + 7 


284 


44.6 


88.8 



38 



Table 11: The Energies (in MeV) and Widths (in keV) of the Electromagnetic 
Mi-Transitions in the (be) Family. 



transition 
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Table 12: The Widths (in keV) of the Radiative Hadronic Transitions in the (bc) 
Family. 



transition 
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Table 13: The Leptonic B c Meson Constant, Calculated in the Different Potential 
Models (the accuracy ~ 15%), in MeV. 



model 
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a 
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13. 44 | 




45 1 
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Table 14: The Leptonic B c Constant, Calculated in the QCD Sum Rules, * is the 
Scaling Relation, in MeV. 
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